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Surface plasmons with ultrasmall optical mode volume and strong near field enhancement can be used to 
realize nanoscale light-matter interaction. Combining surface plasmons with the quantum system provides 
the possibility of nanoscale realization of important quantum optical phenomena, including the 
electromagnetically induced transparency (EIT), which has many applications in nonlinear quantum optics 
and quantum information processing. Here, using a custom-designed resonant plasmon nanocavity, we 
demonstrate polarized position-dependent linewidth-controllable EIT spectra at the nanoscale. We 
analytically obtain the double coherent population trapping conditions in a double- A quantum system with 
crossing damping, which give two transparent points in the EIT spectra. The linewidths of the three peaks 
are extremely sensitive to the level spacing of the excited states, the Rabi frequencies and detunings of pump 
fields, and the Purcell factors. In particular the linewidth of the central peak is exceptionally narrow. The 
hybrid system may have potential applications in ultra-compact plasmon- quantum devices. 



Surface plasmons in metallic nanostructures, which originate from the collective oscillations of free electrons 
in metal 1 , have attracted great interest in light-atom interaction 2 . Compared with other photonic structures, 
plasmon nanostructures provide large and anisotropic Purcell factors caused by an ultrasmall volume of 
optical modes 3 " 7 . Together with the strong near field enhancement near the resonant plasmonic structure, the 
superiority of its Purcell factor has been applied in the nanoscale realization of quantum phenomena, such as 
enhancement and quenching of fluorescence of molecules and quantum dots 8 " 13 , single surface plasmon 
sources 14 " 16 , resonance fluorescence of single emitters 1718 , modification of spontaneous emission 19 , etc. 

Coherent population trapping (CPT) in the stable state of light-atom interaction is the origin of many quantum 
coherence effects. Two kinds of mechanisms for CPT are widely investigated. One is the ordinary two-photon 
resonance leading to the conventional electromagnetically induced transparency (EIT) 20,21 , where the populations 
are coherently trapped in the two lower levels of the three-level A or four-level double-A system 22,23 . Due to 
elimination of absorption and deceleration of the group velocity at the transparent point, EIT has many sub- 
stantial applications in nonlinear quantum optics, such as large Kerr nonlinearity 24 " 26 , four- and six- wave mix- 
ing 27-33 , cross phase modulation 24 ' 26,34 ' 35 , etc. EIT has also played an important role in quantum information 
processing, such as quantum memories 36 " 40 , quantum state transfer 41 ' 42 , quantum gates 16 ' 24 ' 43 , quantum entangle- 
ment 44 ' 45 , etc. 

Another mechanism for CPT is related to the crossing damping between two closely lying upper levels, which 
originates from the interaction with the common vacuum of the electromagnetic field. By involving the crossing 
damping into the trapping condition, the populations can be coherently trapped in the upper levels of a three-level 
V-system or four-level system containing V-configuration, which leads to several spontaneously generated 
coherence (SGC) effects, i.e., spontaneous emission cancelation 46 , quantum beats of population oscillations 47 , 
and narrowing of spectral lines near the trapping condition 48,49 . Crossing damping terms are closely related to the 
Purcell factors in the electromagnetic environment 50 . Plasmonic structures with ultrasmall optical mode volume 
have the advantage of providing large and anisotropic Purcell factors 3 " 7 , thus they can be used to affect the SGC 
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effects of quantum systems. However, the SGC-induced transpar- 
ency in the four-level double-A system is scarcely investigated 51,52 . 

Here, using the custom-designed resonant plasmon nanocavity 
with both large anisotropic Purcell factors and strong near-field excita- 
tion, we demonstrate linewidth-controllable and position-dependent 
polarized EIT spectra of a quantum system at the nanoscale (Fig. 1). 
we first analytically derive the double trapping of coherent populations 
in the double-A quantum system with the crossing damping. One is 
the two-photon resonance with the populations trapped in the two 
lower levels and the other is assisted by the crossing damping, where 
the populations can be coherently pulled into the two upper levels, 
yielding two transparent points in the EIT spectra. Using the evolution 
of the master equation, three-peaked and linewidth-controllable EIT 
spectra are obtained, in agreement with the outcomes from a dressed 
state analysis. The linewidths of three peaks can be effectively modu- 
lated by the level spacing of the upper levels, the Rabi frequencies and 
detunings of pump fields, and the Purcell factors. In particular, the 
linewidth of the central peak is extremely narrow, much less than the 
natural linewidth. This ultracompact plasmon-quantum combined 
system may have potential applications in the nonlinear quantum 
optics and quantum information processing. 

Results 

Double trapping of coherent populations. Consider a closed 
double-A system with two closely lying upper levels |a x ) and \a 2 ) 
and two ground levels \bi) and \b 2 ) (Fig. 1). Both transitions from 
upper levels \ai) with energy fiw ai and \a 2 ) with energy fico a2 to the 
lower level \bi) with energy fico^ are driven by a strong field with 
the frequency v x and detunings defined as An = co fll — co^ — Vi, 
A21 =u>a 2 —cobi — Vi- Simultaneously, a weak probe field with the 
frequency v 2 acts on both transitions from |a x ) and \a 2 ) to \b 2 ) 
with energy fico^ and detunings are A 12 = co ai — co& 2 — v 2> 
A22 = w ai — (Dh 2 — v 2 , respectively. These transitions are also 
coupled by the vacuum modes of electromagnetic fields. In the 
weak coupling region, with the Weisskopf-Wigner approximation, 
the spontaneous decay rate from the upper level |« f ) to the ground 
level I bj) is y^ and the crossing damping between two close upper 
states is k = k x + k 2 . With the dipole and rotation- wave 
approximations, the interacting system can be described by the 
Hamiltonian in the interaction picture 

H fBt =Mii|fli)(fli|+M 2 l|fl2>(fl2|+^(Aii-Ai2)|fe2>(&2| 

-{mu\a l ){b l \+m 2l \a 2 ){b l \+m l2 \a l ){b 2 \ (1) 
+ m 22 \a 2 ){b 2 \+H.c.) 

where = j (2h), for i,j = 1, 2, is the Rabi frequency of the 
coherent field with the complex amplitude 8j coupling the transition 



from the upper level | a t ) to the lower level | bj) with the corresponding 
electric dipole moment %y In the following, the coherent field Sj can 
be provided at the nanoscale by a resonant plasmon cavity. 

We can gain some physical insights underlying the transparency 
windows associated with the phenomena of coherent population 
trapping by finding the eigenstates of the state vectors which are 
decoupled from the vacuum modes to eliminate the fluorescense 
and thus the absorption. The trapping eigenstates corresponding 
to a zero eigenvalue will be the steady solution of equation (7) (see 
Methods), the existence of which requires 

|Q 11 Q 22 -Q 12 Q 21 | 2 -(A 11 -A 12 )(A 11 |Q 21 | 2 + A 21 |Q 11 | 2 ) 

+ f (An - A 12 ) pl±^ |Q 21 1 2 + *»±M |Q n | 2 (2) 

-^atAi-^^i)=o. 

This equation gives i) the conventional two -photon resonance con- 
dition (An = A 12 ), which makes the population trapping in ground 
levels \bi) and \b 2 ); and ii) the SGC-induced trapping condition 

Aii|Q 2 i| 2 + A 21 |Qn| 2 = 0 

and 



— ^11^21 ~ — -Q11Q21 ~~ 0- 

The SGC-trapping leads to the population trapping in levels \ai), 
|a 2 ), and \bi). According to the symmetry of the system, analogous 
trapping conditions for levels \ai), \a 2 ), and \b 2 ) exist. To obtain 
complete SGC-induced transparency, dipole moments between the 
closely lying transitions have to be parallel to generate the maximum 
destructive interference even with anisotropic Purcell factors. When 
double trapping occurs at the same point with A n = A 12 = 
^i2|^ii| 2 /(|^n| 2 + |^2i | 2 )> where (D U = eo fll — co a2 is the upper spa- 
cing, the zero eigenvalue of equation (7) has twofold degeneracy and 
all the superpositions of the two corresponding eigenstates are trap- 
ping states. As shown in the methods, plasmon nanostructures can be 
used to control the decay rates y^ and k { for i, j = 1 or 2 through 
nanoscale anisotropic Purcell factors. 

Predictions of dressed state analysis. To explain the characteristics 
of the numerical absorption spectral profile in the following (Figs. 2, 
3, and 5), now we give the dressed state analysis 53 . Writing the 




Figure 1 | The plasmon-quantum combined system. The schematics of a resonant silver nanocavity, a double-A quantum system with crossing damping 
between two closely lying upper levels, and its dressed states. 
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Figure 2 | Double trapping three-peaked EIT spectra, (a) Absorption and 
(b) dispersion as a function of the probe detuning A 12 in a double- A system 
with parallel dipoles in an isotropic vacuum for upper spacing co 12 = 0.0y 0 
(black), <x>i2 = 2.0y 0 (magenta), <jl> 12 = 4.0y 0 (green), and co 12 = 6.0y 0 
(blue). Q n = Q 21 = 3.0y 0 , Q 12 = Q 22 = 0.03y 0 , = Y zz = 1.0y 0 , and the 
coupling detuning A n = 0.0co 12 . 

interaction Hamiltonian involving the strongly coupled levels \di), 
|a 2 ), and \bi) as 



H = hA u \a 1 )(a 1 \+tiA 21 \a 2 )(a 2 \ 

+ {m u \a l )(b l \+m 2l \a 2 )(b l \+Hx.). 



(4) 



We obtain expressions for its eigenvectors corresponding to 
eigenvalues A l5 /L 2 , and 2 3 
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where Z f = [U, - AnJU- " A 21 ) 2 + (A, - A n ) 2 |Q 21 | 2 + {X { - A 21 ) 2 
| Qi 1 1 2 ] 1/2 . After considering the dissipation and the coupling of the weak 
probe field, we get the equations of motion for the density matrix of 
dressed states and especially the diagonal elements (for more details see 
Supplementary Equation 2). The absorption lines of the probe field are 
related to the fluorescence of spontaneously emitted photons for the 
same transition. Thus for simplicity we only pay attention to the 
modified decay rates of the diagonal matrix elements p D . Di , since they 



determine separately the linewidths of the peaks in the absorption 
spectra and the spontaneously emission spectra, given by the expression 

T D = j I [(y n +y l2 )p n \ 2 (X i -A 21 ) 2 + (y 21 +y 22 )\il 21 \ 2 (X i -A n ) 2 



+ (kj + K 2 )Q.* n Q 21 (It -A n ) - A 21 ) 
+ (k* + k* 2 ) Q„ Q*, {Xi - A u ) {Xi - A 21 )] , 



(6) 



where i = 1, 2, 3. According to the dressed state analysis, the positions 
and linewidths of the absorption lines can be predicted. The linewidths 
of the three-peaked absorption spectra can be effectively controlled by 
the level spacing, the Rabi frequencies and detunings of pump fields, 
and the Purcell factors. 

Linewidth-controllable three-peaked EIT spectra. To fully investi- 
gate the double trapping on varies parameters, we numerically 
calculated the absorption spectra from the ground level \b 2 ) to 
upper levels |a x ) and \a 2 ) through the standard approach based on 
the master equation for the density matrix elements (for more details 
see Supplementary Equation 1). The response of the quantum system 
to the probe field is governed by its polarization P = 
8 o(£2X + £ 2 X*) / 2 and the susceptibility #(v 2 ) is proportional to 

— (^l 2 P ai b 2 ~^ ^22 Pa 2 b 2 )' Dipole moments from the upper levels 

to one ground level in the double-A atomic system are assumed to 
be parallel to demonstrate double-trapping EIT and all dipole 
moments are set to be equal. Dipoles from the upper levels to 
different ground levels are orthogonal, one pair of which are along 
the z axis and the other along the x axis. 

Full numerical calculations for the complex susceptibility as a 
function of the detuning of the probe field are presented in Fig. 2. 
The Rabi frequencies are Q n = Q 21 = 3.0y 0 and Q 12 = Q 22 = 0.03y 0 > 
the detuning of the coupling field is fixed at An = 0.0cl>i 2 , and the 
Purcell factors are set to be Y xx = Y zz = 1.0y 0 - The absorption spectra 
of the imaginary part of the susceptibility (Fig. 2a) show two trans- 
parency channels in an isotropic vacuum with the presence of the 
maximal SGC. The first type of transparency occurring at A 12 = 
0.0co 12 corresponds to the conventional EIT in a double-A system 
under the two -photon resonance condition, which is retained as the 
upper spacing co u varies. The corresponding dark state is the super- 
position of \bi) and \b 2 ), the populations of which are inversely 
proportional to the square of Rabi frequencies. According to the 
population distribution for cl> 12 = 0.0y 0 > the only transparency win- 
dow in the spectrum should be classified as the first type. The second 
one at A 12 = 0.5co 12 , induced by SGC 51 , fulfils the condition in 
equation (3). In contrast to the EIT dark state with empty upper 
levels, the trapping state induced by the destructive interference has 

1 |Oi 2 | 2 :(l + 



the distribution of p flifli : p fl2fl2 : p hihi = 1 + 



in* 



|0 1: 



|^22 1 2 - Thus, the population in \b 2 ) increases with w u . 

1^22 1 

The population distributions from the analytical eigenstates for the 
zero eigenvalue of equation (7) are confirmed by the numerical results. 

The real part of the susceptibility is proportional to the probe 
refraction and its dispersion is related to the group velocity. The 
double-A system provides a positive dispersion as in the usual EIT 
system and leads to subluminal light propagation with V g < c. In 
Fig. 2b, the SGC-induced transparency is accompanied by a remark- 
ably large slope for the real part of the susceptibility and therefore 
slow light propagation compared to the conventional EIT. The group 
velocity at the conventional EIT is invariant but the group velocity at 
the SGC-induced transparency drops as the upper spacing decreases, 
which confirms the different origin of the two types of transparency. 
Similar SGC effects for the group velocity in a V-type atomic system 
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Figure 3 | The linewidth control of three peaks. Absorption as a function of the probe detuning A 12 in a double- A system with varying (a) coupling 
detunings A n , (b) coupling Rabi frequencies Q n , (c) isotropic Purcell factors = T zzy and (d) anisotropic Purcell factors T zz . Other parameters remain 
the same as in Fig. 2 and co 12 = 4.0y 0 . 



have been reported 54 . Removing SGC doesn't change the group velo- 
city at the conventional EIT. 

Since the ground level | bi) is coupled to two closely lying upper 
levels with a strong field, three peaks appear in the absorption spec- 
trum, which characterizes the energy structure of the dressed state 
regardless of SGC. As shown in Fig. 2, the three peaks shift with the 
upper spacing co 12 in the absorption spectra and their positions 
exactly correspond to the eigenvalues of equation (4). Linewidths 
of the central peaks for co 12 = 2.0y 0 , 4.0y 0 , 6.0y 0 are increased as 
0.010y 0 > 0.101y 0 , 0.257y 0 , respectively, in agreement with the predic- 
tion of the dressed state analysis. In this case, the right peak widens 
but the left peak narrows with increasing upper level spacing. The 
destructive interference between parallel dipoles makes the sideband 
linewidths larger and the linewidth of central peak smaller. 

The positions and linewidths of the absorption peaks also depend 
on the coupling detuning A n , the coupling Rabi frequency Qn, and 
the anisotropic Purcell factors as shown in Fig. 3. As the coupling 
detuning grows, the conventional EIT window shifts along the same 
direction due to the two-photon resonance and the peak positions 
move with it, but the SGC-induced transparency point is fixed at the 
probe detuning of A 12 = 2.0y 0 . When A n = 2.0y 0 , the SGC-induced 
trapping condition for the coupling transition is satisfied and the 
transparency window extends within the whole range of the probe 
detuning as shown in Fig. 3a. Deviating from this value of the coup- 
ling detuning in both directions broadens the central peak, in agree- 
ment with a previous report 52 . Furthermore, as the coupling Rabi 
frequencies are increased, the separation between the peaks 
increases, but the two transparent points are fixed since both trap- 
ping conditions are not destroyed (Fig. 3b). Simultaneously, the 
central peak approaches the SGC-induced transparent point and 
becomes narrower and lower. 



Figures 3c, d present the dependence of spectral linewidths on the 
Purcell factors. The Purcell factors don't affect the transparent points 
and the trapping population distributions. We can alter the decay 
rates and the cross damping of all the channels through an isotropic 
vacuum or change either pair of decay rates and cross damping to the 
same ground state by an anisotropic vacuum. When isotropic Purcell 
factors increase, the linewidths of three absorption peaks grow pro- 
portionally and the heights of peaks drop inversely. When only the 
Purcell factor along the z axis T zz increases, similar effects happen 
with reduced magnitude. All numerical calculations of positions and 
linewidths of absorption peaks are consistent with the predictions of 
the dressed state analysis. The above results demonstrate the pos- 
sibility to control the linewidths of the double-trapping EIT spectra 
by a custom-designed resonant plasmon nanocavity, which we now 
present. 

Polarized EIT spectra in a resonant plasmon nanocavity. To realize 
the mechanism for double trapping of populations and allow experi- 
mental investigation of the three-peaked linewidth-controllable EIT 
spectra at the nanoscale, we propose a custom-designed hybrid 
system of the quantum system and the resonant plasmon 
nanocavity. The designed silver nanocavity has several features. In 
order to ensure effective near field excitation of the quantum system, 
the resonance wavelength of the plasmon nanocavity must match the 
wavelength of the pumping transition channels, i.e., the transitions 
from the upper levels |a x ) and \a 2 ) to the lower level \bi). The 
proposed back to back "E"-shaped design in Fig. 4a guarantees 
uniform distributions of electric fields and Purcell factors within 
an area of at least 100 X 100 nm 2 , where the quantum system can 
be located. Furthermore, in this area the near field is strongly 
enhanced, and the absolute values of the Purcell factors are 
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Figure 4 | The resonant silver nanocavity with the near field and Purcell factor distributions, (a) Schematic of a resonant silver nanocavity and (b) its 
absorption. Near field distributions (c) E x and (d) E zy and anisotropic Purcell factors (e) T^ly^ and (f) T zz /y 0 on the xy-plane 70 nm from the metallic 
surface and at the wavelength of X = 600 nm. 



reasonable and their anisotropy is very large, i.e., generally the value 
of T zz /y 0 is much larger than that of T xx /y 0 . Crucially, using present 
nanofabrication techniques, this silver nanostructure can be 
successfully fabricated in the laboratory. 

To demonstrate the above properties, we use Green's tensor 
method 55,56 (see Methods) with the mesh of 10 nm, to characterize 
the resonant silver subwavelength cavity of dimensions x X y X z = 
230 X 200 X 30 nm 3 , as shown in Fig. 4a. Figure 4b displays its 
absorption spectrum showing several peaks, but the largest res- 
onance is at the wavelength of X = 600 nm which matches the 
transitions of |a x ) <-> 1^) and \a 2 ) <-> \bi) of the quantum system. 
It is noted that, for matching different quantum system, the resonant 
wavelength of plasmon nanocavity can be designed on demand. 
Then, we explored the near field distributions and Purcell factor 
distributions at X = 600 nm. It is found that from z = 60 nm to z 
= 140 nm there are strong near field enhancement and large aniso- 
tropic Purcell factors. Increasing the distance from the metallic sur- 
face to the quantum system decrease the near fields and the Purcell 
factors approach 1.0. Most importantly, in the xy-plane, there is an 
area of about 100 X 100 nm 2 with almost uniform electric fields and 
Purcell factors as shown in Figs. 4c-f. Generally, the larger the electric 
field E x and the Purcell factor T xx /y 0 , the smaller E z and T zz ly Qy and 
vice versa. The anisotropy of electric fields guarantees both the polar- 
ized excitation of the quantum system and linewidth control using 
the near field (or Rabi frequencies) in EIT spectra, and the anisotropy 
of Purcell factors provides the effective control of linewidths of three 
EIT peaks. 



In the xy-plane of z = 100 nm, we choose three points 0, 1, and 2, 
at which the quantum system can be located. They are shown in 
Fig. 4a, point 0 is at the center and points 1 and 2 are set as front 
and right 50 nm from the center. We let two sets of dipoles be 
orthogonal, one is along the x-axis and the other z-axis. First, we 
let the dipoles of |a x ) <-> | bi) and \a 2 ) <-> \ b\) be excited by the optical 
near field E x of the silver nanocavity and the dipoles of |a x ) <-> \b 2 ) 
and \a 2 ) <-> \b 2 ) are scanned by the probe field. It is found that the 
positions and linewidths of the three peaks in EIT spectra are not very 
sensitive to the location of the quantum system, shown as the green, 
blue and red curves in Fig. 5a (the right part of which shows the 
linewidth details of the central peak), which comes from the uniform 
design of electric field and Purcell factor distributions. Next, we let 
the channels of excitation and scanning be exchanged. At points 0 
and 1, the three-peaked EIT spectra are not presented due to very 
weak near fields (not shown in the figure). However, for point 2, 
when the quantum system is excited by the near field of the z dir- 
ection, there is an obvious change in positions and linewidths of three 
peaks, shown as black curves in Fig. 5a. Thus for the near field 
excitation of different directions, polarized EIT spectra are clearly 
observed. 

Finally, we calculate the EIT spectra of the quantum system at the 
center for the different distance from the metallic surface at z = 
80 nm (point 3), 100 nm, 120 nm (point 4), and 140 nm (point 5), 
respectively. The results, shown in Fig. 5b with the right part of 
amplified central peaks, indicate that the peaks and widths of three 
peaks are very sensitive to the distance of the quantum system from 
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Figure 5 | Polarized EIT spectra in a resonant plasmon nanocavity. EIT 

spectra of the quantum system as a function of the probe detuning A 12 for 
(a) different locations on the xy-plane of z = 100 nm and (b) different 
distances for z = 80, 100, 120, 140 nm. Right insets are the amplified center 
peaks. 

the metallic surface due to the strong modification of near fields and 
Purcell factors. To sum up, position-dependent linewidth-control- 
lable polarized EIT spectra of the designed quantum system is rea- 
lized at the nanoscale in a plasmon nanocavity. 

Discussion 

It is noted that, in the above SGC-induced trapping condition [equa- 
tion (3)] of a quantum system with crossing damping, it is required 
that the two sets of dipoles are strictly parallel. This is difficult to fulfil 
in a general atomic system, thus experimental realizations of this 
kind of system are seldom investigated 52 ' 57 . One possible scheme to 
observe the phenomena experimentally is using the Rubidium 85 
hyperfine structure 52 . The two ground states \bi) and \b 2 ) are 5Si /2 , 
F = 3 and 5S 1/2 , F = 2. The upper state 5P 1/2 , F = 3 is coupled to 
5D 3/2 , F = 4 by a strong field. Thus the upper state splits into two 
dressed states corresponding to \a^) and |a 2 ), according to the 
dressed state analysis. Then the atomic system can fulfill the trapping 
condition that dipoles between \a^) <-» \b t ) and \a 2 ) <-> \b t ) (for i = 1, 
2) are parallel. Another potential quantum system to form the above 
four-level system is dual CdSe/ZnS/CdSe nanocrystals 58,59 . Their 
photoluminescence peak positions lie in the optical window, such 
as X = 600 nm, and their lifetime is on the order of nanoseconds. The 
position of each state in both the nanostructures isolated by the wide 
band gap barrier can be independently tuned, which leads to a design 
of the desired four-level system by changing materials and sizes. 
Furthermore, altering the thickness of the barrier layer can control 
the coupling between the two upper states. 

In summary, we have theoretically investigated the EIT of the 
four-level double-A quantum system and its nanoscale realization 
in the resonant plasmon nanocavity. We have first analytically 
derived the double trapping conditions of populations, leading to 
two transparent points in the EIT spectra. By evaluating the master 
equations, we have numerically obtained the three-peaked line- 
width- controllable EIT spectra, and found the positions and line- 
widths of peaks in agreement with the predictions of the dressed 



state analysis. Further numerical investigations indicate that the line- 
widths of these peaks can be well modulated by the level spacing of 
two upper levels, pump Rabi frequencies and detunings, and iso- 
tropic and anisotropic Purcell factors. Finally, the combined system, 
composed of the custom-designed resonant plasmon nanocavity and 
the quantum system, has been proposed to realize the three-peaked 
linewidth-controllable polarized EIT spectra at the nanoscale. This 
study bridges the fields of traditional quantum optics and plasmonics 
associated with the strong near field enhancement and ultrasmall 
optical mode area and will benefit low photon number quantum 
nonlinear optics and quantum information processing. 



Methods 

The equations of motion for the state vectors. With the Weisskopf-Wigner 
approximation, the dissipation is incorporated via the equations of motion for the 
state vectors, A^t), A 2 (t), B^t), and B 2 (t): 



dt 



M(t)- 



2 

y 2 i+y 2 2 N 



A 2 (t)+in ll B l (t)+in l2 B 2 (t), 



A 2 (t)- 



^Ai(f)+in 2 iBi(f)+/n 22 B 2 (0, 



(7) 



B l (t) = iQ* l A l {t) + iQ* 2l A 2 (t), 

B 2 (t) = i{ An - An )B 2 (t) + n* l2 A x (t) + iQ* 22 A 2 (t) . 



The upper spacing co 12 is much less than the dipole transition frequencies, 
guaranteeing that the dipole moments between the upper levels and one of the ground 
levels interact with the common electromagnetic vacuum. This leads to crossing 
damping terms k = K\ + k 2 and k* = k\ + k 2 , which represent the quantum 
interference between the spontaneous emission pathways and strongly depend on the 
properties of the vacuum. In a general vacuum with anisotropic Purcell factors, the 
intersection angle between the dipole moment fx^ and the x axis is Oy, for i, j = 1, 2. If 
we let the y axis be the quantum axis, the damping terms are given by 

Vij = ( r ** cos2 % + r ^ sin2 %) and Kj = (F xx cos 0 y cos 6 2j + 

l^lll I Mil I 

T zz sin By sin 6 2 j), for i,j = 1,2, where T zz ly 0 and T xx ly 0 are the Purcell factors in the* 



and z directions, and y 0 



is the decay rate in a vacuum for the 



transition between |a x ) and \bi). According to these expressions, large values of 
crossing damping, which lead to destructive interference, appear in a vacuum only 
with nearly parallel dipoles 48 or in an anisotropic vacuum even with orthogonal 

dipoles 19 ' 49 . 

Green's tensor method. In the design of resonant plasmon nano structure, Green's 
tensor method for solving the near field and anisotropic Purcell factors is used 55,56 . 
When an arbitrary shaped subwavelength structure with the electric permittivity 
6(r, co) embedded in the homogeneous bulk material with eo (co), the field E(r) at any 
point r satisfies the Lippmann-Sch winger equation 



E(v)=E\v)+k 2 f dr'G°(r, r', co)e s (r', co)-£(r'), 
Jv 



(8) 



where V denotes the clusters subspace, e s (r, co) =6(r, co) — Go(co), and Green's tensor 
in a three dimensional system is 



G°(r, r', Q)=[I 



1 - ik 0 R 



RR J^^ (9) 



where R = |R| = |r — r'| and k 2 Q = /c 2 eo(co). The needed anisotropic Purcell 
factors 1718 ' 60 can be expressed as the Green's tensor coefficients 
^zz/Yo = 3 A fll &1 ImG(r) zz and T xx /y 0 = 3X aibl lmG(r) xx in the x and z directions, 
where G(r) = G° (r) + k 2 \ v dr'G 0 (r, r', co)e s (r , co) -G(r ) . 
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